Thermal rectification effects of multiple semiconductor quantum dot junctions 
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Based on the multiple energy level Anderson model, this study theoretically examines the thermo- 
electric effects of semiconductor quantum dots (QDs) in the nonlinear response regime. The charge 
and heat currents in the sequential tunneling process are calculated by using the Keldysh Green's 
function technique. Results show that the thermal rectification effect can be observed in a multiple 
QD junction system, whereas the tunneling rate, size fluctuation, and location distribution of QD 
significantly influence the rectification efficiency. 



I. INTRODUCTION 

Due to global warming and environment issues, it is 
important to developing the new energies to reduce the 
total amount of CO2 in the earth's atmosphere. Electri- 
cal power generated by the mechanical force of wind and 
ocean tides is one of several types of green energy. Apart 
from that, electrical power generated from solar energy 
transformation also creates clean energy. Researchers 
have used many technologies to develop materials with 
high values of figure of merit (ZT) for potential appli- 
cations of solid-state thermal devices. These materials 
can be utilized to generate electrical power through the 
temperature bias resulting from solar heating or other 
method heating. [1-9] 

Researcher have proposed several methods for enhanc- 
ing ZT.[2] One of them is to reduce system dimen- 
sions. [3] A zero dimension QD system was predicted to 
be more pronounced for the enhancement of thermoelec- 
tric efficiency in the dimension reduction. [9] Neverthe- 
less, previous studies focus on the linear response regime 
of AT /To <C 1, where AT and To are the temperature 
bias and equilibrium temperature, respectively, of two 
side electrodes. [10] The properties of thermal devices in 
the nonlinear response regime AT /To ~ 1 have impor- 
tant applications such as thermal rectifiers and thermal 
transistors. A thermal rectifier is crucial for the heat cur- 
rent storage. Records of thermal rectification date back 
as early as 1935, when Starr discovered that copper ox- 
ide/copper junctions can exhibit a thermal diode behav- 
ior. [11] The thermal rectification effects have been the- 
oretically predicted to occur in one dimensional phonon 
junction systems. [12-15] Because phonons are heat flow 
carriers in refs [12-15], their applications are restricted 
to the heat manipulation and heat energy storage. Using 
charge carriers (electrons or holes), we can manipulate 
heat currents and charge currents. The Peltier effect with 
the application of coolers makes it possible to manipulate 
the heat current using applied voltage. Its reverse pro- 
cess is so called the Seebeck effect used to create electrical 
power. 

Scheibner and coworkers have, recently, employed a 
single metallic QD junction to measure electrochemi- 
cal potential in the linear response regime. [16] Such a 



small electrochemical potential yielded by temperature 
bias is not sufficient to demonstrate the thermal rec- 
tification effect. [17] Besides, the operation temperature 
of a metallic QD is very low because its charging ener- 
gies are much smaller than the thermal energy (fesT) = 
25 meV). Semiconductor QD junctions with large in- 
terdot Coulomb interactions and energy-level separation 
within each QD allow semiconductor QD thermal recti- 
fiers to operate at much higher temperatures. Previous 
studies of refs[16,17] only consider a single QD junction, 
which is useful for the case of low QD density. To scale 
up QD thermal devices, it is necessary to consider the 
case of high QD density. 

This study demonstrates that a system of multiple 
semiconductor QDs embedded in an amorphous insula- 
tor with low heat conductivity connected to the metallic 
electrodes exhibits the thermal rectification effect in the 
nonlinear response regime. This study also clarify the 
thermal power in the linear and nonlinear regimes. Al- 
though the mechanism of thermal rectification for QD 
junctions is similar to that of charge current, the heat 
current is yielded by temperature bias and electrochem- 
ical potential. The very nonlinear relation between the 
applied temperature bias and the electrochemical poten- 
tial leads to that it is not straightforward to reveal the 
relation between the heat current and the applied tem- 
perature bias. 



II. FORMALISM 

The following multi-level Anderson model [18] 
describes the Hamiltonian of the metal/quantum 
dots/metal double barrier junction in Fig. 1 ; 
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where a k a s (a,k,o,p) creates (destroys) an electron of mo- 
mentum k and spin a with energy in the j3 metal- 
lic electrode. The term d\ a (c^ l<7 ) creates (destroys) 
an electron with the ground-state energy Ei of £th QD 
, while Ug and Ugj describe, respectively, the intradot 
Coulomb interactions and the interdot Coulomb interac- 
tions . The term Vk,p,t describes the coupling between 
the band states of electrodes and the QDs. The Hamilto- 
nian of Eq. (1) assumes that energy level separations be- 
tween the ground state and the first excited state within 
each QD is much larger than intradot Coulomb interac- 
tions Ug and thermal temperature kgT. Therefore, there 
is only one energy level for each QD. This study ignores 
the interdot hopping terms due to the high potential bar- 
rier separating QDs. 

Using the Keldysh-Green's function technique, [19] the 
charge and heat currents leaving electrodes can be ex- 
pressed as 

Je = ^E/ de 7i (e)ImGl a (e)f LR (e), 

Q = yX! / de le (e)ImGl a (t)(t-E F -eAV)f LR (e), 

where /iij(e) = f L (e) - /r(c), the transmis- 
sion factor is 7/ (e) = r^S+r^) ' fL(R)( £ ) = 
l/(exp^~^ L ( R ^^ kBTL ( R ^ + 1) is the Fermi distribution 
functions for the left (right) electrode . The chemi- 
cal potential difference between these two electrodes is 
related to the bias difference hl — = eAV. The 
temperature difference is Tl — Tr = AT, and Ep is 
the Fermi energy of electrodes. T^ ^e) and r^ ^j(e) 
[^ej = 27r^ k \Ve^^\ 2 S(e — e k )] denote the tunneling 
rates from the QDs to the left and right electrodes, re- 
spectively. Notations e and h denote the electron charge 
and Plank's constant, respectively. For simplicity, we 
assumed that these tunneling rates are energy- and bias- 
independent. Equations (2) and (3) have been employed 
to reveal the thermal properties of a single-level QD in 
the Kondo regime. [20,21] Here, this analysis focuses on 
the heat current in the Coulomb blockade regime. Previ- 
ous studies [22,23] provide an expression of the retarded 
Green function 

+ JV '-'E e -B,-^-n„ + j r,' (1) 

where n denotes the number of QDs. The term H m de- 
notes the sum of Coulomb interactions seen by a particle 
in dot I due to other particles in the dot j(j ^ £), which 
can be occupied by zero, one or two particles. The proba- 
bility factor p m denotes the probability of such configura- 
tions. For a three-QD system (£ ^ j ^ j'), there are nine 



(3 x 3) configurations, and the probability factors become 

Pi = a j a j'' P2 = bjClj,, p 3 = djbj,, Pi = Cjdy, p 5 = Cj,CLj, 

Po = bjbji, p'j — Cjbji, ps = Cjibj, and pg — CjCj>, where 
Oj = 1 - (N jt a + N jt _ a ) + cj, bj = {N jt<r + Nj^) - 2c„ 
and Cj = (nj.- a nj. a ) is the intradot two particle corre- 
lation function. Nj i(r is one particle occupation number. 
The interdot Coulomb interaction factors are II i =0, 

n 2 = u tlj , n 3 = u e , f , n 4 = 2U itj , n 5 = 2U e , f , 

n 6 = U e ,j + U e ,j>, n 7 = 2U e ,j + U e ,j>, n 8 = 2U Lr + U ej , 
and Ilg = 2Ugj + 2U(j'. The sum of probability factors 
p m for all configurations equals to 1, indicating the fact 
that G\ CT (e) satisfies the sum rule. The imaginary part of 
each resonant channel is Ye, — (J~t,L + ^e,R))/2 resulting 
from the coupling between the QDs and the electrodes, 
where the real part of self energy is ignored. A self energy 
that ignores the effects of electron Coulomb interactions 
is adequate in the Coulomb blockade regime, but it does 
not capture the Kondo effect. 

According to the expression of the retarded Green's 
function in Eq. (4), we need to know the occupation 
numbers A^ CT (A^_ CT ) and Nij — eg. These numebrs 
can be obtained by solving the following equations self- 
consistently: 

Nt a = _ j dg r^/^ + r^^) 

J it T e , L + r e ,R ■ ' 1 ' 

The values of Ne, a and q are restricted between and 1. 
The expression of two particle retarded Green function 
of Eq. (6) is 

GUe) = N t _„ £ - _ _ — - — 

m—i 



III. RESULTS AND DISCUSSION 

Although the thermal rectification effect of a single 
QD junction with two levels has been theoretically in- 
vestigated in ref.[17], it is difficult to realize the physi- 
cal condition of asymmetrical coupling between the QD 
and the electrodes, where the excited state states have 
symmetrical tunneling rates, but the ground state has 
an asymmetrical tunneling rate. Meanwhile, ref[17] did 
not take into account the crucial intralevel Coulomb in- 
teractions. To determine the direction-dependent heat 
current, we considered two sets of parameters (the cases 
of T L > T R and T R > T L ): (a) is fi L = E F + eAV, 
Hr = E F , T L = T + AT, and T R = T , and (b) 
is /il = Ep, fiR = Ep — eAV ,T L = T and T R = 
T + AT. The average temperature of two side electrodes 
is (T L + T R ) /2 = T + AT/2, which is increasing steadily 
with AT. Such a physical condition for average temper- 
ature is more readily realized than that of reference [17] 
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from experimental point of view. [24] A functional ther- 
mal rectifier requires a good thermal conductance for the 
case of Tl > Tr, but a thermal insulator for the case 
of T R > T L . Based on Eqs. (2) and (3), a QD thermal 
rectifier requires not only highly asymmetrical coupling 
strcngthes between the QDs and the electrodes, but also 
strong electron Coulomb interactions between the dots. 
To investigate the thermal rectification behavior, we nu- 
merically solved Eqs. (2) and (3) for three QD junctions 
for using various system parameters. 

The energy level of dot t is Eg = E F + ottAE, where 
— 1 < on < 1. The term AE = 200r represents the 
energy level fluctuation of QDs. For simplicity, the in- 
tradot Coulomb interactions Ut = 600T are fixed, al- 
though the QD size fluctuation leads to the intradot 
Coulomb interaction variation. According to Eq. (3), the 
direction-dependent heat currents are attributed to Ay 
and Ni i<7 = Nt t _ a = N(. The Seebeck effect for the open 
circuit ( J e = 0) determines the electrochemical potential 
eAV. This electrochemical potential forms in response to 
the current generated by the temperature gradient. This 
electrochemical potential is the Seebeck voltage (Seebeck 
effect). The thermal power is S — AV/ AT. Once AV is 
solved, we then use Eq. (3) to compute the heat current. 

Figure. 2 shows the direction-dependent heat cur- 
rent (Q) and rectification efficiency (t]q) as functions of 
AT for various values of Tar, while keeping T b(c).l = 
Fb(c).r = r. We adopted the interdot Coulomb interac- 
tions U e ,j = 300r and k B T = 25T. The energy levels of 
dots A, B, and C were chosen to be Ea = E F — AE/5, 
E B = Ep + 2AE/5 and E c = E F + AE/5. The heat 
currents have positive and negative signs corresponding 
for the cases of Tl > Tr and Tr > Tl, respectively. 
The former is the heat current from the left electrode 
to the right electrode Qlr, while the later is the heat 
current from the right electrode to the left electrode 
Qrl- The rectification efficiency is r\Q = (Qlr(AT) — 
\Q RL {AT)\)/Q LR {AT). The rectification effects of heat 
currents in Fig. 2(a) are clearly distinguishable from Fig. 
2(b). When dot A is fully blocked (T A r = 0), the QD 
junction exhibits a much higher t\q and the heat currents 
decrease dramatically. This serious suppression implies 
that the heat currents through dots B and C are small 
because their resonant energy levels are far away from 
the Fermi energy level due to a large Uac and U ab ■ The 
results in Fig. 2 indicate that precisely manipulating the 
coupling strengthes between the dots and the electrodes 
is crucial to the optimization of thermal rectifiers. 

Figure. 2 illustrates the homogenous interdot Coulomb 
interactions. However it is difficult to achieve homoge- 
nous interdot Coulomb interactions in the multiple QD 
layout. The interdot Coulomb interaction fluctuation 
mainly arises from the size variation and the inhomoge- 
nous position distribution. Therefore, we examined the 
variation effect of Uac ° n heat currents. Figure. 3 
shows the heat current and rectification efficiency as func- 
tions of AT for various values of Uac at ^ar = and 
Tal = 2r. The other parameters are the same as those in 



Fig. 2. In this case, large interdot Coulomb interactions 
suppress the heat currents. This behavior is similar to 
the Coulomb blockade effect on the charge current. 22,23 
Nevertheless, a strong Uac leads to higher rectification 
efficiency. Therefore, a higher efficient thermal rectifica- 
tion effect appears in the high density QD system. In 
the such a case of strong intradot and interdot Coulomb 
interactions, the charge current exhibits clear plateaus 
with respect to applied voltage. However, the structure 
of heat currents in Fig. 3 does not exhibit these plateaus 
with respect to temperature bias. This is because tem- 
perature bias broadens the electron distribution. [22] 

In Figures (2) and (3) we illustrate the case of Ec < 
Eb- Figure. 4 shows that Ec influences the heat 
current and t]q. The other parameters are the same 
as those in Fig. 3. When Ec changes from Ec — 
E F + AE/5 to E c = E F + 4AE/5, the heat current 
for E c = E F + 2 AE/5 is less than that for E c = 
Ep + 3 AE/5. This indicates that the separation be- 
tween the energy levels of QDs and the Fermi energy 
level is not a unique factor to influence the heat current 
magnitude. When Ec increases to Ec — Ep + AAE/5 
from Ec — Ep + 3AE/5, the heat currents and rectifica- 
tion efficiency change slightly. All of these complicated 
behaviors are the result of the very nonlinear relationship 
between the heat currents and electrochemical potentials, 
which are strongly influenced by the electron Coulomb 
interactions and the energy levels of each QD. 

The following expression of heat current, Q = Q B + 
Qc, further clarifies the results of Fig. 4: 

Q B /(j B n) = (1- N b )*{{1-2Na){E b - E F )f LR (E B ()4) 

+ 2N A (E B + Uab - E F ) * / LH (£ B + Uab) 

Qc/(lcn) = {l-N c )*{{l-2N A ){Ec-E F )f LR {Ec) 

+ 2Na(E c + Uac - E F ) * f LR {E c + U AC )(5) 

Equations. (7) and (8) are the results of setting q 
to zero due to very large intradot Coulomb interactions. 
Further, a delta function replaces the Lorentzian func- 
tion of resonant channels for T fcsTo in Eq. (3). In 
addition, Eqs. (7) and (8) only consider the interdot 
Coulomb interactions Uab and Uac since there is one 
particle occupation in dot A. As for Na, Nb, Nc and 
eAV in Eqs. (7) and (8), we use full numerical solution 
(Eqs (2), (5) and (6)) to solve them. Figure. 5 shows the 
heat current as functions of AT for Ec = Ep + 3AE/5. 
The other parameters are the same as those in Fig. 4. The 
lot lines given in Eqs. (7) and (8) match the full solutions 
very well. Qlr and Qrl are the nonlinear functions of 
AT. Figure. 5(b) shows the Qb and Qc calculated by 
results of Eqs. (7) and (8) for Qlr and Qrl- Based on 
Fig. 5(b), Qc is the primary determinant of the heat 
currents of Fig. 5(a). The factors of Na and Jlr in- 
fluence the behavior of Qlr with respect to AT. Only 
the factor of Jlr influences the behavior of Qrl- This 
is because of AT-independent Na = 0.45 in the case of 
Qrl- So far, we have analyzed the heat current and t\q. 
From experimental point of view, it is easy to measure 
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the electrochemical potential yielded by the temperature 
bias although it is not straightforwardly related to the 
heat current in this system. 

Figure. 6 shows the electrochemical potential of differ- 
ent values of UbTq. The other parameters are the same 
as those in Fig. 5. The dot and dashed lines represent 
the linear functions of AT, whereas the solid lines are 
nonlinear with respect to AT. Although the asymmetri- 
cal electrochemical potentials (Fig. (6a) and (6b)) reveal 
some information about thermal rectification, it is dif- 
ficult to determine the rectification efficiency from the 
measurement of electrochemical potentials. This is true 
for both linear and nonlinear responses due to the very 
nonlinear relation between Q and AV. Therefore, it is 
difficult to judge the thermal rectification effect from the 
experimental results done in ref[16]. The large electro- 
chemical potential created by the temperature bias (AT) 
is useful for generating electrical powers. The results of 
Fig. 6 also imply that the charge current rectification 
with respect to temperature appears in the closed circuit 
(J e 7^ 0). An electrochemical potential is an essential 
quantity for charge carriers to deliver heat current, it is 
different from the phonon carriers in a phonon junction 
system. 12-15 To further illustrate the strong nonlinear 
relationship between electrochemical potential and heat 
current, we show the heat current and rectification ef- 
ficiency in Fig. 7, where the curves exhibit one to one 
correspondence with those in Fig. 6. We note that tjq 
is enhanced when temperature fc_eT decreases. A com- 
parison of Fig. 6 and Fig. 7 demonstrates the strong 
nonlinear relationship between Q and AV. 

Finally, figure 8 calculates the thermal power as a func- 
tion of fceTo for different values of AT. The other pa- 
rameters are the same as those in Fig. 5. As expected, 
the temperature bias has a significant influence on ther- 



mal powers at low fcsTo, but not at high temperatures. 
Note that the thermal power values increase with increas- 
ing temperature bias. Thermal power plays a significant 
role in determining the figure of merit. [10] Some thermal 
devices such as solid state coolers are in the nonlinear re- 
sponse regime. 2 Therefore, the previous studies of linear 
response regime underestimate the thermal power values 
when AT/To > 1. [10,25] 



IV. CONCLUSIONS 

This theoretical study reports the thermal rectifica- 
tion effects (TRE) of multiple QDs embedded into an 
amorphous insulator (phonon glass) with low heat con- 
ductivity. Results show that the asymmetrical tunneling 
rates and strong interdot Coulomb interactions signifi- 
cantly influence the TRE. Recent research indicates that 
it is possible to precisely manipulate the size and po- 
sition of multiple germanium QDs in S1O2 amorphous 
insulator. [26] This experiment indicates that it is possi- 
ble to realize the system shown in Fig. 1 to examine 
the thermal rectification at high temperatures. To avoid 
the phonon heat current, which can seriously suppress 
the heat rectification effect arising from electrons, scan- 
tunneling- microscopes (STMs) can replace one of the 
electrodes since the vacuum layer between the STM and 
amorphous insulator can blockade the phonon heat cur- 
rent. [27] 
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Figure Captions 

Fig. 1. Schematic representation of multiple semicon- 
ductor quantum dots embedded into an amorphous in- 
sulator with low heat conductivity sandwiched between 
metallic electrodes. 

Fig. 2. Heat current and rectification efficiency as 
functions of AT for the tunneling rate variation of dot 
A. The heat currents are in the units of Qo = T 2 /(2h). 
Such units are used through out this article. 



Fig. 3. Heat current and rectification efficiency as 
functions of AT for different interdot Coulomb interac- 
tions of Uac at Far, = and Tal = 2I\ The other 
parameters are the same as those in Fig. 2. 

Fig. 4. Heat current and rectification efficiency as 
functions of AT for different energy levels of Ec and 
Uac — 300r. The other parameters are the same as 
those in Fig. 3. 

Fig. 5. Heat current as functions of AT for Ec = 
Ep + 3AE/5. The Other parameters are the same as 
those in Fig. 4. 

Fig. 6. Electrochemical potential as functions of AT 
for different values of fcsT . The other parameters are 
the same as those in Fig. 5. 

Fig. 7. Heat current and rectification efficiency as 
functions of AT for different values of UbTq. The curves 
are one to one correspondence to those of Fig. 6 

Fig. 8. Thermal power as functions of fee To for the 
different values of AT. The other parameters are the 
same as those in Fig. 5. 
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